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TRANSONIC ROTOR TIP DESIGN USING NUMERICAL OPTIMIZATION 


Michael E. Tauber and Ronald G. Langhi* 
Ames Research Center 


SUMMARY 


The aerodynamic design procedure for a new blade tip suitable for operation at 
transonic tip speeds is illustrated. For the first time, three-dimensional numeri- 
cal optimization was applied to rotor tip design, using the recent derivative of the 
R0T22 code, program R220PT. Program R220PT utilized an efficient quasi-Newton 
optimization algorithm. Multiple design objectives were specified. The delocaliza- 
tion of the shock wave was to be eliminated in forward flight for an advance ratio 
of 0.41 and a tip Mach number of 0.92 at ip = 90®. Simultaneously, it was sought to 
reduce torque requirements while maintaining effective restoring pitching moments. 
Only the outer 10/t of the blade span was modified; the blade area was not to be 
reduced by more than 'i%. The goal was to combine the advantages of both swept-back 
and swept-forward blade tips (configurations which had been studied previously). A 
planform that featured inboard sweepback was combined with a swept-forward tip using 
a taper ratio of 0.5. Initially, the R0T22 code was used to find by trial and error 
a planform geometry which met the design goals. This configuration had an inboard 
section with a leading edge sweep of 20® and a tip section swept forward at 25®; in 
addition, the airfoils were modified. With this tip shape as a starting point, 
program R220PT was used to improve the planform and airfoil section design. Signif- 
icant improvements, which met all design objectives, were achieved in 2.5 CPU hours 
on a CRAY X-MP computer. The optimized tip has an inboard sweepback of 19.7°, a tip 
swept forward at 34.1®, 50% of taper, and modified airfoil sections. 


INTRODUCTION 


The inevitable increase in forward flight speeds of conventional helicopters 
has resulted in higher rotor tip Mach numbers. However, forward flight speed can be 
limited by the increased wave drag caused by the presence of large regions of super- 
sonic flow on the outer part of the advancing blade. Another undesirable feature 
resulting from this supersonic flow can be a large increase in the noise emanating 
from the rotor blade. Strong shock waves forming on the blade tip can propagate 
efficiently into the far field, producing intense low-frequency noise. To lessen 
these drag and noise increases, rotor tips designed to reduce shock strengths have 
been extensively studied (ref. 1-6). (In fact, most new high-speed helicopters have 
swept-back rotor tips.) 
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Although swept-back rotor tips have been successfully used in recent years, 
such a tip configuration is generally not optimum. Since the tip is normal to the 
oncoming free stream in the second quadrant, the torque requirements can exceed 
those for a rectangular blade in that position (refs. 6 and 7). Reference 6 theo- 
retically shows that a swept-forward tip offers both torque and acoustic advantages 
over a swept back one. However, swept-forward lifting surfaces have a well-known 
tendency toward angle-of-attack divergence instability which can produce severe 
vibration (although it may be possible to restrain the divergence by using aniso- 
tropic composite materials). To avoid the structural and dynamic problems of the 
swept-forward tip, it was decided to Investigate tip shapes which would be aerody- 
namically stable, yet combine the best features of both forward and aft sweep. The 
theoretical study was made using the R0T22 transonic, full-potential, three- 
dimensional rotor code and its recent derivative, R220PT, which is suitable for 
optimization. 

The next section briefly discusses the development of R0T22 and R220PT, 
including the main features and the simplifying assumptions of both codes. A 
description of the application of the codes to the design of a new, multiply-swept 
blade tip follows. The modification of both the tip planform and the airfoil 
sections in the tip region will be demonstrated on a hypothetical modern lifting 
rotor blade with an initially rectangular planform and a supercritical airfoil. 


BRIEF DESCRIPTION OF ANALYSIS 


Program R0T22 

The earliest codes developed for calculating the high-speed flow over a rotor 
blade solved the transonic small disturbance equation (refs. 8-10). The R0T22 code 
represented the first solution of the full potential equation for the flow over a 
rotary wing. Since R0T22 has been previously documented (refs. 11 and 12), its main 
features and assumptions will be only briefly reviewed here. 

Beyond the assumptions inherent in the potential equation, the simplifications 
embodied in the model are: a single-blade analysis, a simplified mathematical 

treatment of the vortex sheet behind the blade, an approximate wake-induced inflow, 
and a quasi-steady flow field. Because the full -potential equation and the exact 
tangency boundary condition are used, subsonic and transonic flow over thick air- 
foils with blunt leading edges and at angle of attack can be modeled realisti- 
cally. There are no inherent geometric limitations; airfoil section and chord can 
be varied, together with pitch angle, angle of attack, twist, and flapping angle. 

The code is inviscid and contains a relatively simple, curved, kinematically deter- 
mined near-wake. Wake-induced inflow can be approximately modeled by judicious 
placement of vortices in the flow field. Provisions exist for two straight vor- 
tices. Thus, the tip vortex shed from a preceding blade can be approximately simu- 
lated. The strength of the vortices is determined from the circulation on the blade 
and the vortex core size is based on the correlation of reference 13. 
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The governing equations (potential and energy) are rendered quasi-steady by 
setting all the time derivatives of the perturbation potential equal to zero; how- 
ever, the effects of rotation are maintained with the quasi-steady assumption. 
Although the time history is removed from the solution for forward flight, the equa- 
tions and the boundary conditions retain a timelike dependence by being a function 
of rotor position. The assumption of quasi-steady flow permits computing the flow 
field about the blade at one azimuth angle. This speeds the computation by nearly 
two orders of magnitude in comparison with a fully unsteady calculation. (In an 
unsteady computation, the azimuth angle can be increased only a fraction of a degree 
between time steps as the blade advances into the airstream (ref. 14).) The effect 
of the quasi-steady assumption was assessed in reference 15. Here, computed surface 
pressure distributions around the entire rotor disk were compared with measured 
values near the tip on both straight- and swept-tip model rotors to advance ratios 
of 0.45 and to a maximum tip Mach number of 0.91. The agreement between calculated 
and measured pressures was found to be good. In fact, at the important 90° blade 
azimuth position, the theoretical and experimental values compared very well to an 
advance ratio of 0.5, corresponding to a tip Mach number of 0.94. 

While many sources of helicopter noise have been identified (ref. 16), impul- 
sive noise was found to dominate the noise spectrum at high speed (ref. 17). In 
reference 17, Schmitz and Yu explained the cause of high-speed impulsive noise as 
the propagation of shock waves occurring on the advancing blade into the acoustic 
far-field, a mechanism which they labeled "delocalization." Thus, the accurate 
prediction of velocities in the vicinity of and beyond the blade tip is critical to 
the determination of transonic noise propagation into the acoustic far-field. It 
had been shown previously (eq. (1), ref. 6) that rotor-generated noise caused by the 
formation of shock waves on the blade propagated most effectively when the blade was 
at the 90° azimuth position because the sonic surface was closest to the blade 
tip. The unperturbed sonic surface is located at 


S = 



2 2 
V cos 



- u sin i|» 


( 1 ) 


where S is measured from the center of rotation, M^.^^ is the rotational Mach 
number, y is the advance ratio, and i|» is the azimuth angle. R0T22-computed and 
experimentally measured velocities were extensively compared in this part of the 
flow field (ref. 18). In the experiment, laser velocimeter measurements were made 
of the three orthogonal velocity components in the neighborhood of the tip of a 
nonlifting, hovering model rotor at tip Mach numbers from 0.85 to 0.95. The mea- 
surements were made in three vertical planes, one located two-thirds of a blade 
chord inboard of the tip, one at the tip, and one at the unperturbed sonic circle. 
Each extended one chord above the blade. The R0T22-calculated velocities agreed 
very well with the LV measurements, as did the prediction of the strength and posi- 
tion of the shock waves. Since the angle of incidence of the advancing blade tip in 
the vicinity of 90° is usually small, the R0T22 code can be used with the reasonable 
confidence to predict the important tip-region flow field. The same utility extends 
to R220PT which is derived from R0T22. 
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Program R220PT 


Linking R0T22 with minimization algorithm QNMDIF has resulted in program 
R220PT, which retains all of the capabilities of R0T22 but has been drastically 
reorganized to meet the requirements of numerical optimization. These requirements 
include evaluating multiple design points in the same run (typically the fully 
advancing and fully retreating blade positions), matching prescribed pressure or 
Mach number distributions at arbitrary computational span stations, perturbing the 
blade sections or planform (or both), and computing the flow field for each pertur- 
bation as efficiently as possible. 


Quasi-Newton Optimization Method 

QNMDIF is a general-purpose, unconstrained, quasi-Newton minimization program 
using finite difference approximations to first derivatives. The quasi-Newton 
algorithm seeks the minimum of a quadratic model to an objective function of N 
variables by stepping in a descent direction derived from gradient information and 
from a matrix that approximates second derivatives. The second derivative informa- 
tion improves as optimization proceeds. If the objective were truly quadratic, the 
quasi-Newton method theoretically would locate the minimum in N steps. Any 
desired constraints must be included in the objective as penalty functions or be 
imposed implicitly by the construction of the problem. 

Typically, the accuracy of the flow calculations is far less than the available 
machine precision. The current implementation of QNMDIF is effective even with 
objective functions having less than full machine precision. It uses a scheme to 
estimate appropriate finite differencing intervals, which contributes to successful 
optimization of nonanalytic functions such as those being evaluated (refs. 19 
and 20). 


Design by Optimization 

The basic concept of aerodynamic design by optimization is well documented 
(refs. 20 and 21); specify target aerodynamic characteristics and seek to match 
them by reshaping the lifting surface geometry. The difference between the computed 
flow quantities and the collection of targets constitutes the objective function to 
be minimized. Each evaluation of this function involves at least one flow field 
computation, which is itself an iterative process. Reliable optimization needs a 
minimum of four or five significant digits in the objective. R220PT monitors the 
precision of the objective by evaluating and printing it at every iteration of each 
flow solution. The precision required by the optimizer typically demands a higher 
level of flow convergence than is normally considered adequate for R0T22. A maximum 
correction criterion on the order of 1.E-07 is usually needed for optimization. 

Objective function evaluations dominate the cost of a run, and the number of 
evaluations required is directly proportional to the number of design variables 
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used. R220PT accepts as few as one input station per planform discontinuity in 
defining the blade geometry. Since section modifications are applied at the input 
stations, the sparse blade definition leads to an economy of design variables. 

High convergence criteria and multiple function evaluations of this type place 
a heavy load on a computer's central processing unit (CPU). Even with the well- 
vectorized flow solver common to both R220PT and R0T22, one flow solution from start 
using a 120 X 16 X 48 mesh takes about 250 CPU seconds for 400 iterations on a CRAY 
X-MP. Both R220PT and QNMDIF provide for starting from existing flow solutions and 
optimization data, respectively. This allows a run to be interrupted and resumed 
without a loss of information. The flow solutions within an optimization run use a 
neighboring flow field as a starting guess, with a limit of 200 iterations for any 
one solution imposed after the first function evaluation. Many of the evaluations 
require far fewer iterations. 


BASIC CONFIGURATION AND DESIRED MODIFICATIONS 


This study postulates a hypothetical modern rotor blade which is representative 
of those used on high-speed helicopters. The blade was assumed to have the follow- 
ing characteristics (ref. 6): an aspect ratio of about 10, a linear twist variation 

from hub to tip of -10°, and an NLR-1 supercritical airfoil having a maximum 8.71t 
thickness-to-chord ratio. A rotational tip Mach number of 0.65 (725 ft/sec) and an 
advance ratio of 0.41 were used; the corresponding forward flight speed was 
90.6 m/sec (176 knots). The maximum tip Mach number, at 4» = 90°, was about 0.92. 
The cyclic pitch angle, 0^^, of the advancing blade was described by 

9 = 10.7 - 11.2 sin i> (2) 

c 

and was -0.5° at i(i = 90° and 1° at both ^ = 60° and 120°. The conditions were 
considered to be representative of typical high-speed flight. This geometrically 
simple blade with a rectangular tip exhibited strong delocalization at these flight 
conditions. 

Only the outer W% of the blade (0.977 chord length) was modified (ref. 6). 

The amount of tip taper was constrained so that the decrease in total blade area was 
less than 3%- For structural reasons, it was desirable to decrease the airfoil 
thickness-to-chord ratio no more than 20% so that the value at the blade tip would 
not be less than 7%. The following performance improvements over a rectangular 
blade were stipulated for the azimuth angle range from 60° to 120°; (1) delocaliza- 

tion was to be eliminated, (2) the torque requirements were to be reduced, (3) the 
restoring pitching moments were to be at least as large, and (4) the thrust coeffi- 
cients were not to be decreased significantly. The R0T22 code was used to find, by 
trial-and-error , an improved tip configuration. Once the feasibility of meeting the 
design objectives had been proven, R220PT was used to achieve additional performance 
improvements. 
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Planform Modifications 


Two types of tip-planform modifications were investigated. These consisted of 
taper and various combinations of fore and aft sweep. Initially, the tip taper was 
based on the results of a previous design study (ref. 6) which had established the 
effectiveness of using 50% taper on the outer 10% of the blade in reducing or elimi- 
nating delocalization. The resulting 2.8% blade area reduction was within the 
desired 3% limit. When the R220PT code was used, tip taper was made a design varia- 
ble; however, the total area reduction was constrained to 3%. The sweep modifica- 
tions were applied after dividing the outer 10% of the blade into two equal segments 
of about one-half chord length each. The two segments were then swept at different 
angles. Sweeping the inboard half of the tip forward and the outboard half aft 
moved the aerodynamic center of the tip forward of that of the remainder of the 
blade, and produced the same undesirable pitching moment behavior characteristic of 
the entirely swept forward tip (ref. 6). It was decided, therefore, to concentrate 
on using sweepback inboard and to sweep the tip forward. The leading edge sweep 
angles which were used in ROT22 computations ranged from 10® to 25® for the inboard 
part, and from -15% to -25% for the outboard part. When R220PT was used, the lead- 
ing edge sweep angles of both portions of the tip were made design variables. 


Airfoil Modifications 

The contours of the NLR-1 airfoil were modified at selected spanwise blade 
stations on the outer 10% of the blade by adding geometric shape functions to the 
original ordinates of the airfoil. The shape functions were added linearly to the 
baseline profile's upper and lower surface to achieve the desired aerodynamic 
improvements. The surface perturbations were faired linearly to zero at the 90% 
spanwise blade station. 

In the R0T22 code, three types of shape functions were used to modify the 
airfoil sections in the tip region. These consisted of an exponential, a sine, and 
a droop function (ref. 22). Combinations of these three functions had previously 
been used successfully to design proposed two-dimensional (ref. 21) and three- 
dimensional (ref. 23) rotorcraft lifting surfaces. The exponential function predom- 
inantly affects the leading-edge region and is useful either in achieving efficient 
supercritical flow or in modifying leading-edge pressure gradients to improve stall 
progression. The sine function is especially useful for weakening shock waves by 
changing local curvature of the airfoil at selected chordwise locations; it can be 
used more than once at each section. The droop function, which was added symmetri- 
cally to the upper and lower surfaces, produces a vertical displacement of the 
forward region of the profile; it is useful in controlling leading-edge pressure 
gradients and in eliminating shocks on the lower surface in the leading-edge region. 

In the R220PT code, airfoil section modifications were made by adding Wagner 
functions (ref. 24) to the original airfoil's ordinates. The Wagner functions were 
more practicable when optimization was used since fewer variables were needed to 
describe the airfoil modifications than if shape functions had been used. In 
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practice, however, some flexibility was lost. Since the number of design variables 
had to be limited to control computation time, the Wagner functions were only 
applied to the upper surface of the blade. Had shape functions been used, for exam- 
ple, the upper and lower leading edge region could have been modified simultaneously 
with the droop function. 


RESULTS AND DISCUSSION 


To illustrate how the codes were used to design a blade tip for more effective 
and quieter operation at transonic tip speeds, the flow was calculated on the rec- 
tangular and modified tip configurations. Since the blades all developed lift, a 
longitudinal vortex was placed lO^t inboard of the tip and one chord below the blade 
to simulate a vortex shed from a preceding blade. Using a mesh having 120 cells 
chordwise, 48 spanwise, and 16 vertically resulted in 33 computational planes posi- 
tioned normal to the blade surface. Nonlinear spanwise grid spacing was used to 
improve resolution in the tip region; seven of the computational stations were 

placed on the outer 10/t of the blade, with another seven beyond the blade tip. The 

outermost grid plane was 0.1556 blade radii (about 1.5 chords) beyond the tip. The 

combination of the 0.41 advance ratio and the rotational tip Mach number of 0.65 
resulted in the idealized sonic surface being 0.1285 blade radii from the tip at 
4» = 90® (from eq. (1)). The sonic surface was, therefore, well within the 
computational domain at the 90® blade azimuth position where delocalization was most 
likely to occur. 

The flow field, and especially delocalization, can be graphically presented by 
plotting contours of constant Mach numbers from the solutions. Up to 20 contour 
lines can be shown at chosen Mach numbers, both on the blade surface and off the 
surface in the plane of the blade. Since any desired contour can be highlighted, 
the use of a heavy dashed line for the sonic line and dotted ones for the remaining 
contours permitted easy identification of the extent of delocalization. 

First, the flow field on the blade with the previously listed double-cranked 
tips was calculated using the R0T22 code. In this manner, a tip geometry which met 
the design objectives was identified. This tip geometry was then used as a starting 
point for optimization to illustrate how the R220PT code can be used to improve the 
initial design. 


Program R0T22 Tip Design 

Initially, the flow field about the basic unmodified rectangular blade was 
calculated using the ROT22 code. At the 90® azimuth position, the cyclic pitch was 
-0.5®; the advance ratio of 0.4l and the rotational tip Mach number of 0.65 resulted 
in a tip Mach number of about 0.92. Although all calculations were made using the 
entire blade, for clarity only enlarged views of the tip region flow field will be 
shown. Contours of constant Mach number on the upper (higher-speed) surface of the 
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blade and in the plane of the blade are shown in figure 1 . The Mach number contours 
generally began at 0.8 and went to the peak value in 0.05 intervals, although the 
more critical transonic values were plotted in increments of 0.02. 

Impulsive noise occurs when the shock wave terminating the supersonic flow 
concentrated in the tip region of the advancing blade extends into the far-field, 
beyond the blade tip. In the blade-attached coordinate system of the present calcu- 
lation method, delocalization would correspond to an open sonic line. This delocal- 
ization is clearly illustrated in figure 1 . The heavy dashed contour represents the 
sonic line. The large size of the supersonic region beyond the tip, which even at 
its narrowest point is nearly half a blade chord wide, is apparent. The shock 
forming on the blade propagated into the flow field beyond the tip, perturbing the 
sonic surface significantly from the location predicted by equation (1). (Note also 
that about one chord inboard of the tip on the blade surface there are two super- 
sonic zones, separated by a narrow strip of subsonic flow. This flow pattern is a 

result of the presence of a double shock, an undesirable phenomenon occurring on 

supercritical airfoils operating at off-design conditions.) 

The inviscid coefficients of thrust, C.p, torque, Cq, and pitching moment, C|^^, 
were calculated for the entire blade at advance angles from 60° to 120°, using the 
cyclic pitch variation of equation (2). The force and moment coefficients for the 

rectangular blade will be used as reference values for comparison with those for the 

blades having modified tips. The results for the rectangular blade are presented in 
table 1 . 

The R0T22 code was used to calculate the flow on the entire blade at the 90° 
azimuth position where delocalization was most likely to occur. The tip taper was 
fixed at 50 % while the sweep of the two portions of the tip section were systemati- 
cally varied. The values ranged from 10° to 25° sweepback for the inboard part and 
-15° to -25° of (forward) sweep for the tip. Approximately five or six combinations 
of sweep were investigated in increments of 5°. Each computer run required about 
one minute of CPU time on the CRAY X-MP. As the sweep angle of the two sections was 
increased to 25° and -25° (the largest values considered for each part) delocaliza- 
tion was eliminated. However, a small pocket of high Mach number flow formed aft of 
the "V" where the two swept sections joined. Since the flow was accelerated as the 
sweep angles increased, it was decided to reduce the sweepback of the inboard por- 
tion to 20°. Although delocalization reoccurred, the supersonic region beyond the 
tip was only 0.1 of a chord length wide at its most narrow part (fig. 2). However, 
the peak Mach number on the blade, aft of the V, was still between 1.40 and 1.45, 
resulting in a very strong shock followed by a weaker one (fig. 3). Although the 
strong shock was limited to a very small region on the blade, the resultant wave 
drag raised the torque to that of the rectangular blade. The thrust and pitching 
moments were nearly the same for both blade configurations. It was therefore, 
necessary to modify the NLR-1 airfoil sections on the entire tip to try to eliminate 
delocalization. In addition, the shock strengths had to be reduced to decrease the 
wave drag responsible for the high torque. 

The airfoil modifications were made by applying the previously described expo- 
nential, sine, and droop shape functions. The functions were used at each of the 
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seven stations with the sine function being used twice at all except the tip. The 
initial choice of coefficients (ref. 23) for the functions was based on previous 
experience with modifying the NLR-1 airfoil section (ref. 6). Approximately 15 com- 
puter runs were made during which the most important of the 74 shape function coef- 
ficients were manually varied. The resulting modifications eliminated delocaliza- 
tion (fig. 4) and reduced the torque by over 20?. The thrust and pitching moments 
were only minimally affected. Not surprisingly, it was necessary to make the air- 
foil about 10? thinner at the V; this reduced the peak Mach number by almost 
0.15. The corresponding shock strengths have also been greatly reduced on both the 
upper and lower surfaces (fig. 5). 

Additional flow calculations were also made over the azimuth arc from 60* to 
120®. The present preliminary configuration exhibited no delocalization, required 
10? to 20? less torque, and matched the thrust and pitching moment behavior of the 
rectangular blade. The Mach number contours and pressure distributions for 
ij) = 60° and it» = 120® are shown in figures 6 and 7, and 8 and 9 , respectively. 

The peak Mach numbers on the blade were 1.35 at 4» = 60° and 1.30 at 4» = 120®. 
Next, the application of optimization will be illustrated in an effort to improve 
the design further. 


Program R220PT Tip Design 

Numerical optimization was applied to the 20®/-25® planform of the final hand- 
modified configuration described above, with the original NLR-1 airfoil section 
retained. All optimization computations were made at the 90® blade azimuth posi- 
tion. R220PT was run with 48 computational span stations, such that the V created 
by the intersection of the cranked portions of the blade was aligned between two 
stations. This placement was chosen because the steep spanwise velocity gradients 
existing there cause the flow field calculations to oscillate if a computational 
station is at or very near the V. Blade section modifications were effected by use 
of Wagner shape functions, combinations of which resemble airfoil surfaces 
(refs. 20, 22, and 24). 

Optimization began with 15 variables, consisting of the two tip planform sweep 
angles and the tip chord, and two sets of the first six Wagner functions. The 
Wagner functions were applied to the airfoil upper surface at the V and the blade 
tip; the functions were linearly faired to zero at the 90? blade radial station. 

The object was to reduce the peak Mach numbers in the region from 5? of the span 
inboard of the tip to 6.7? outboard. Constraints were imposed on the total blade 
force coefficients of thrust, torque, and pitching moment via penalty functions to 
achieve the desired aerodynamic characteristics. The tip area was held constant by 
adjusting the chord at the V as the tip chord varied. The initial optimization 
step took 2350 CPU seconds on the CRAY X-MP (running under COS 1.13). The resulting 
Mach number contour plot showed that the sonic zone delocalization had been elimi- 
nated. The crank angles had changed to 19.7° and -34.1®, while the tip chord 
remained essentially unchanged at 50? taper. The force coefficients were accepta- 
ble, and the flow computations converged smoothly. 
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Subsequently, the flow solutions began to oscillate as the forward sweep of the 
tip section Increased further. Because the configuration resulting from the first 
optimization step provided smooth convergence, It was decided to fix the planform at 
the angles of 19.7® and -34.1®, and to concentrate the limited amount of time 
remaining upon Improving the airfoil sections. Eight to twelve Wagner function 
variables were used to modify the airfoil sections. After 23 additional optimiza- 
tion steps, the design objectives had been met. These final 23 steps required 
6800 CPU seconds, bringing the total computing time for the optimization to Just 
over 2.5 CPU hr. (Although the optimization was only done at the 90® blade azimuth 
position, the 2.5 hr of CPU time should not be considered excessive. A previous 
design optimization of a two-dimensional rotorcraft airfoil (ref. 21) required a 
roughly comparable amount of computing time. (Approximately 30 hr of CPU time on a 
CDC 7600 computer were used. However, the CDC 7600 was an order of magnitude slower 
than the CRAY X-MP computer used in the present study.) 

Next the computed performance characteristics of the optimized ( 19.7®/-34. 1®) 
tip will be discussed and compared with the starting tip shape of 20® /-25® and, 
finally, with the rectangular blade. The Mach number contours for the 90® advancing 
blade position are shown in figure 10. Although the peak value is 1.35 (or 0.1 
higher than for the 20®/-25° tip), the width of the subsonic region beyond the tip 
is 50/5 greater for the optimized tip. The optimized tip also had a 15/5 greater 
restoring pitching moment than the 20®/-25® tip. However, the torque was 6/5 higher, 
partly since the maximum surface Mach number was higher and partly because the lower 
leading edge lip shock was still present on the inboard, swept back portion 
(fig. 11). 

The performance of the optimized tip over the remainder of the advancing blade 
arc was also calculated, although it was not included in the optimization process. 
The Mach number contours at = 60® are shown in figure 12. The maximum value is 
1.45, or again 0.1 higher than for the starting tip shape. The higher Mach number 
caused a strong shock to form aft of the V (fig. 13). The strong shock increased 
the torque requirement by 1 % over that for the starting shape; however, the pitching 
moment improved by 19?. At the 120® blade position, the peak Mach number was again 
0.1 higher (see fig. 14) than for the configuration that was used to start the 
optimization. The shock on the inboard, swept back part was very strong (see 
fig. 15). However, the torque was only 2.5? higher for the optimized tip, while the 
pitching moment value improved by 10?. 

By comparing the Mach numbers on the optimized configuration (figs. 10, 12, 
and 14) with those of the starting one (figs. 2, 6, and 8), it can be seen that they 
are consistently higher, especially aft of the V. The higher Mach numbers resulted 
from the greater change in direction that the flow on the optimized tip experienced 
at the V, and is proportional to the included angle of the V. For the optimized 
tip, the included angle was 126®, compared to 135® for the tip shape at the start of 
optimization. Had viscous effects been present in the codes, the optimized config- 
uration might have had a shallower included angle between the two tip portions. 
Otherwise, the strong shocks would likely induce boundary layer separation which 
would have resulted in a large increase in drag and torque. 
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Finally, the torques and pitching moments were compared for the optimized 
configuration (19°/-34.1®) and the start-of-optimization shape (20®/-25°) with two 
previously studied tip shapes (ref. 6). The tip shapes of reference 6 both had the 
same 50% taper; one had 30® leading edge sweepback, while the other had 20° of 
leading edge forward sweep. All results were normalized by those for a rectangular 
blade, therefore including the rectangular tip in the comparison also. However, a 
more important reason for presenting normalized values is to suppress viscous 
effects which are important in determining torques. The normalized torque coeffi- 
cients for blades with the four tip configurations are shown in figure 16 for 
advancing blade angles from 60° to 120°. Note that the 20° swept forward tip had 
the lowest torque and the 30° swept back tip the highest from iji = 70° to 120°. The 
double-cranked tips, however, always required less torque than the rectangular 
blade, with the 20°/-25® being the lower one. The optimized tip's torque was 1% 
below the rectangular blade at i|j = 60°, l6/( below at 90°, and b% less at 120°. It 
was in the restoring pitching moment comparison that the optimized tip fared best, 
however, while the disadvantage of the 20° swept forward tip was most obvious 
(fig. 17). The pitching moments of the blade with the optimum tip were 36/S higher 
than the rectangular configuration at 4) = 60° , 10/S higher at 90°, and 7/S higher at 
120°. These values were from lOlS to 20% higher than the pitching moments of the 
20°/-25° tip; therefore, the design objectives had been met. The blade area reduc- 
tion of 2.8% was within the 3/S limit that was specified. The thinnest airfoil 
section was at the tip (see fig. 15) and had a maximum thickness-to-chord ratio of 
7.04/S, corresponding to a decrease of 19?. (Since the thrust coefficients of the 
modified blade differed from those of the rectangular blade by only a few percent 
over the advancing blade arc from 60° to 120°, they were not shown.) 


CONCLUDING REMARKS 


Three-dimensional numerical optimization was used for the first time in aerody- 
namic design of a high-speed rotor tip. A tip planform geometry consisting of a 
swept-back section and a forward-swept tip, combined with 50? taper, meets all the 
design objectives. Twenty-four computer optimization cycles were required using 
from 8 to 15 design variables. The 2.5 CPU hr of computing time on a CRAY X-MP was 
considered moderate for a three-dimensional aerodynamic optimization calculation. 

When comparing the improvements in the design achieved through optimization, 
several limitations should be kept in mind. Since numerical optimization in three- 
dimensional flow is not a fully automated process, the dominant limitation here was 
time, both in man-hours and in active CPU time. For instance, time constraints 
limited the application of optimization to the 90° blade azimuth position only. Had 
other blade positions been included in the study, such as 60°, 120°, and 270°, the 
final configuration could have been altered substantially. Also, had more time been 
available, the Wagner functions could have been applied to the lower surfaces of the 
airfoil sections, instead of only the upper surfaces. Thus, the wave drag and 
resultant torque could have been further reduced. Therefore, the optimized config- 
uration should not necessarily be considered an optimum design, although all the 
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design objectives have been met. Rather, the results illustrated here should be 
considered an example of the successful application of the R220PT code to the compu- 
terized aerodynamic design process. 
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TABLE 1.- FORCE AND MOMENT COEFFICIENTS FOR 
UNMODIFIED RECTANGULAR BLADE 


— 

deg 

9c 

deg 

Cip 

Cq 

^M 

60 

1.0 

5.251x10"3 

9.280x10“5 

-1.504x10*^ 

75 

-.1 

4.677x10"3 

5.547x10*5 

-2.084x10“^ 

90 

-.5 

4.328x10"3 

3.751x10*5 

- 2 . 307 x 10 *^ 

105 

-.1 

4.833x10’3 

4.758x10"5 

- 2 . 128 x 10 ’J 


1.0 

5.629x10'3 

8.106x10*5 

- 1 . 677 x 10 ”^ 
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Figure 1.- Mach contours on rectangular tip with NLR-1 airfoil section, i|i = 
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Figure 2.- Mach contours on 20®/-25® swept tip with NLR-1 airfoil section, if 



Figure 3.- Pressures on unmodified NLR-1 airfoil sections; sweep 20®/-25°, 
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Figure 3-- Concluded. 
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Figure 4.- Mach contours on 20®/-25“ swept tip with "hand-modified” airfoil 

sections, i|» = 90®. 
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Figure 5.- Pressures on "hand-modified" airfoil sections; sweep 20°/-25®, = 90°. 
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Figure 5.- Concluded. 
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Figure 6.- Mach contours on 20* /-25® swept tip with '• hand-modified" airfoil 

sections, il) = 60®. 
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Figure 7.- Concluded. 
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Figure 8.- Mach contours on 20® /-25® swept tip with "hand-modified" airfoil 

sections, <l> = 120°. 




Figure 9.- 


Pressures on "hand-modified" airfoil sections; 


sweep 20“/-25®, = 120”. 
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Figure 10.- Mach contours on optimized tip, sweep 19.7®/-34. 1® , i|) = 90 



Figure 11.- Pressures on optimized tip; sweep 19.7®/-34. 1® , = 90®. 
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Figure 11.- Concluded. 
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Figure 13.- Pressures on optimized tip; sweep 19.7®/-34. 1®, = 60® 
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Figure 15.- Pressures on optimized tip; sweep 19.7°/-34. 1® , 






36 











CQ/CQ RECTANGULAR 


Figure 16 


O UNMODIFIED SECTIONS, PLANFORM 30 deg SWEPT BACK 

O HAND-MODIFIED SECTIONS, PLANFORM ANGLES 
20/-25 deg 

A OPTIMIZED SECTIONS, PLANFORM ANGLES 19.7/-34.1 deg 

+ UNMODIFIED SECTIONS, PLANFORM 20 deg SWEPT 
FORWARD 



Normalized 


torque comparison for blades with 


four modified tips. 
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□ UNMODIFIED SECTIONS, PLANFORM 30 deg SWEPT BACK 
O HAND-MODIFIED SECTIONS, PLANFORM ANGLES 
20/-25 deg 

A OPTIMIZED SECTIONS, PLANFORM ANGLES 19.7/-34.1 deg 
-I- UNMODIFIED SECTIONS, PLANFORM 20 deg SWEPT 
FORWARD 



Figure 17.- Normalized pitching moment comparison for blades with four modified 

tips. 
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